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Abstract 

Rayleigh- Taylor instability in a compressible fluid is reconsidered. The density is allowed 
to vary with pressure under the barotropy assumption. For the case with equal speeds of 
sound in the two superposed fluids, in order to give a non-trivial compressibility correction to 
the Rayleigh- Taylor growth rate, the compressibility correction is calculated to 0{g'^/k'^a^). 
To this order, compressibility effects are found to reduce the growth rate. 
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1 Introduction 



The instability of the interface between two fluids having different densities and accelerated 
towards each other is called the Rayleigh- Taylor instability (Rayleigh [T], Taylor [2J). For 
the case of two superposed uniform fluids separated by a plane interface, the amplitude of a 
small disturbance which is periodic in the horizontal interfacial plane will vary in time like 
e""*, where 

n = ■ gk 1 

P2 + Pi 



while pi and p2 are, respectively, the densities of the lower and the upper fluids, k the 
horizontal wavenumber of the disturbance, and g the gravitational acceleration. (1) shows 
that if P2 > Pi, the interface is unstable. Chandrasekhar [3] and Hide included the effect 
of viscosity while Bellman and Pennington [5] showed that the effect of surface tension T at 
the interface in the linear problem is to produce a critical wavenumber kc, given by 

1 

2 

(2) 



g{p2- Pi) 



so the interface is unstable or stable according to whether the wavenumber k is less than 
or greater than kc- Rayleigh- Taylor instability plays a crucial role in inertial conflnement 
fusion (Petrasso [6]) and gravitational fusion in stars (Burrows [7]). 

The effect of compressibility on Rayleigh- Taylor instability was considered by Vander- 
woort [8] and Plesset and Hsieh |9] but the results were in disagreement with each other 
(Shivamoggi [TQj). Several attempts have since been made to clarify the role of compressibil- 
ity (Plesset and Prosperetti [11], Bernstein and Book [12]), but the picture is not completely 
clear. In view of the inconsistency, as noted by Plesset and Prosperetti [TT], in assuming 
the constancy of both density and speed of sound a by Vanderwoort [S], in this paper, we 
reformulate Vanderwoort's [S] development by relaxing the constant density assumption and 
allow density to vary with pressure under the barotropy assumption. The resulting disper- 
sion relation turns out to be different from the one given by Plesset and Prosperetti [llj 
in general, but agrees with the latter to 0{g/ka'^). However, the latter compressibility cor- 
rection disappears when the speeds of sound in the two fluids become equal. In order to 
get a non-trivial result for the latter case, in this paper we proceed further and give the 
compressibility correction to next order, namely 0{g'^ /k'^a^). 



2 Governing Equations 

Consider the static state of inflnite fluid arranged in horizontal layers in a uniform gravi- 
tational fleld y, the density p and the pressure p are functions of the vertical coordinate z 
satisfying the equation of hydrostatic equilibrium - 

dp 

-T- = -P9- (3) 

az 
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To determine the character of this equihbrium, consider the evolution of small disturbances 
5p in p, 5p in p and v in velocity on this static equilibrium state. The linearized equations 
governing these disturbances are - 



continuity: 



motion: 



d5p 
~dt 



+ V • Vp + pV • V = o 



dv 

= -V5p - g6p U 



(4) 



(5) 



and following Vanderwoort [8j, we use that barotropic condition (this is what indeed Van- 
derwoort [8] used without explicitly mentioning it) relating p and p - 



d6p ^ f d6p 



+ V ■ Vp 



(6) 



a being the speed of sound. 

Assuming the disturbances to vary with x, y and t according to Q^ikix+k2y)+nt obtain 
(Vanderwoort [8J) - 



P 



(P + nVa2) 



Dw 



n^pw + gk'^ 



D + 



where, 



a"^/ \ (fc^ + 11? /a?) 



d 



w = o 



(7) 



V =< u,v,w > and D = — . 

dz 

The relevant boundary conditions at an interface supporting discontinuities in p and a 
are the kinematic jump condition - 



Aw = o 

and the dynamic jump condition that follows from equation (7) (Vanderwoort [8]) 



P 



Dw 



+ gk'^A 



P 



(k^ + nya^) 



{w)o = o. 



(9) 



3 The Rayleigh- Taylor Configuration 

Consider now the case of two semi-infinite fluids separated by a horizontal plane interface at 
z = o with the pressure gradient in each fluid given by equation (3). Following Vanderwoort 
[8], we assume the speed of sound in the two fluids to be constant so that p is a linear 
function of p. However, the density under the barotropy assumption coupled with equation 
(3) in each fluid is now given by (Mathews and Blumenthal [13] ) 



-Dp= — -. 



(10) 



p 

It may be noted that the constant density assumption coupled with the constant sound 
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speed assumption used by Vanderwoort [8], as noted by Plesset and Prosperetti [TT], is 
inconsistent. 

Using equation (lo), equation (7) becomes 

D-w - (I) Dn^-{k- + '^y = o (11) 

which agrees with the one given by Plesset and Prosperetti |llj . 
One may seek a solution to equation (11) of the form 



which then leads to 



from which. 



w ~ e^" (12) 



Assuming the disturbances to be bounded at infinity and invoking the jump condition 
(8), we obtain 



W = < aoz ^ (15) 



Using (15), the jump condition (9), gives the dispersion relation 



n 



P2 



k"^ + m? /al 



q2 - 



Pi 



+ 



P2 



Pi 



k'^ + 'n?/a\ k'^ + n^/a\ 



o. 



(16) 



Here, subscripts 1 and 2 refer to the upper and the lower fluids, respectively. 
Introducing 



/c' 1 + 71^ /k'^a? 



(16) may be rewritten as 



n 



gk{a2 - ai] 



ai/ii - a2fi2 

It may be noted that the dispersion relation given by Plesset and Prosperetti [11 
not agree with (18) - in the present notation, the one given by [H] is 



n 



gk{p2 - Pi) 



Pi/Pi - P2/P2 

which is different from (18). However, (19) turns out to agree with (18) to O (g/ka^) . 



(17) 

(18) 
does 

(19) 
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4 Small Compressibility Limit 

The complexity of (16) (or (18)) makes physical interpretations harder without a numerical 
calculation. It is useful to consider (16) in the small compressibility limit - {g/ka?) ^ 1. 
In this limit, on noting from (17) that 



± 



+ 



+ 1 



(20) 



or 



± 1 + 



n 



o 



+ 



(21) 



where hq is the growth rate in the incompressible limit {{g/ka^) 0) 



no 



Igk 



P2- Pi 
P2 + Pi 



we obtain from (18), 



n 



1 + 



pip2g 



Hpi + P2y \aj 



k^af. 



(22) 



(23) 



which agrees with the one given by Plesset and Prosperetti |llj . 

In the special case ai = a2, (23) implies that compressibility effects disappear from 
this problem to 0{g/ka?). We now need to proceed further and include in (23) terms to 
0(//A;V). 

First, on noting from (17) that 



2ka{. 



± 



1 + 



+ 



2i2 



2k^al2 Sk^af 



+ 



k^al 



we obtain from (18), with Oi = 02 = a, 



n 



P1P2 



+ 



(24) 



(25) 



(25) shows that, in the special case ai = a2, compressibility effects have a stabilizing character 
on this problem - in this special case, the general belief that fluid compressibility provides a 
sink for the potential energy released at the interface appears to be valid. 



5 Discussion 

In this paper, we have reconsidered Ray leigh- Taylor instability in a compressible fluid. The 
density is allowed to vary with pressure under the barotropy assumption. For the case with 
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equal speeds of soTind in the two superposed fluids, in order to give a non-trivial compress- 
ibility correction to the Rayleigh- Taylor growth rate, we have calculated the compressibility 
correction to O {g"^ / k'^a^) . To this order, compressiblity effects are found to reduce the growth 
rate. In this special case of equal speeds of sound, the fluid compressibility appears to provide 
a sink for the potential energy released at the interface. 
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